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Antiproton-induced reactions on nuclei at the beam energies from hundreds MeV up to several GeV provide 
an excellent opportunity to study interactions between the antiproton and secondary particles (mesons, baryons 
and antibaryons) with nucleons. The antiproton projectile is unique in the sense that most of the annihilation 
particles are relatively slow in the target nucleus frame. Hence, the prehadronic effects do not much influence 
their interactions with the nucleons of the nuclear residue. Moreover, the particles with momenta less than 
about 1 GeV/c are sensitive to nuclear mean field potentials. This paper discusses the microscopic transport 
calculations of the antiproton-nucleus reactions and is focused on three related problems: (i) antiproton poten- 
tial determination, (ii) possible formation of strongly bound antiproton-nucleus systems, and (ili) strangeness 


production. 
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I. MOTIVATION 


It is difficult to produce antiproton beams. However, 
antiproton-nucleus interactions have attracted experimental- 
ists and theorists since about 30 years when the KEK and 
LEAR data appeared. Since this time significant progress has 
been done to describe these data on the basis of optical and 
cascade models. Still, antiproton interactions inside nuclei 
need to be better understood. One example is the antiproton- 
nucleus optical potential. According to the low-density theo- 
rem, it can be expressed as 


2 
Vopn = FE fin. (1) 


where at threshold Vs = 2my, E> =~ my, fpp =~ (—0.9 + i0.9) 
fm [1]. Being extrapolated to the normal nuclear density pọ = 
0.16 fm’, Eq.(1) predicts the repulsive antiproton-nucleus 
potential, ReVop ~ 75 MeV. In contrast, the p-atomic X-ray 
and radiochemical data analysis [2] favors the strongly attrac- 
tive antiproton-nucleus potential, ReVopt ~ —100 MeV in the 
nuclear center. Thus the pA optical potential is not a simple 
superposition of vacuum PN interactions. The strongly at- 
tractive pA potential is consistent with the Relativistic Mean 
Field (RMF) models and has a consequence that a nucleus 
may collectively respond to the presence of an implanted an- 
tiproton. The formation of strongly bound p-nuclei becomes 
possible [3, 4]. 

Another very interesting aspect is p-annihilation in the nu- 
clear interior. This results in a large energy deposition > 2my 
in the form of mesons, mostly pions, in a volume of hadron- 
ic size ~ 1 — 2 fm [4, 5]. After the passage of annihilation 
hadrons through the nuclear medium a highly excited nuclear 
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residue can be formed and even experience explosive multi- 
fragment breakup [5, 6]. The annihilation of an antiproton 
at Pab <~ SGeV/c on a nuclear target gives an excellent 
opportunity to study the interactions of secondary particles 
(pions [7], kaons and hyperons [8], charmonia [9, 10]) with 
nucleons. This is because most annihilation hadrons are slow 
(y < 2) and have short formation lengths. Thus their interac- 
tions are governed by usual hadronic cross sections. 

Over the last decades, several microscopic transport mod- 
els have been developed to describe particle production in 
pA interactions [6, 7, 11-13]. Nowadays there is a renais- 
sance in this field, since the antiproton-nucleus reactions at 
Pla ~ 1.5 — 15 GeV/c will be a part of the PANDA experi- 
ment at FAIR. The most recent calculations are done within 
the Giessen Boltzmann-Uehling-Uhlenbeck (GiBUU) model 
[14-16] and within the Lanzhou quantum molecular dynam- 
ics (LQMD) model [17, 18]. In the present paper, I will report 
some results of the GiBUU calculations for p-nucleus inter- 
actions at Piab = 1.5 — 15 GeV/c. 


Il. GIBUU MODEL 


The GiBUU model [19, 20] solves a coupled set of kinetic 
equations for baryons, antibaryons, and mesons. In a RMF 
mode, this set can be written as [21, 22] 


(py! pO, aps (ea; 4 mð m D= = fix, p *) 


=G. (2) 


where œ = 1,2,3, w = 0,1,2,3, x = (tr) j = 
N, Ñ, 4, 4, Y,¥, n, K, K etc.. f;(x,p*) is the distribution 
function of the particles of sort j normalized such that the 
total number of particles of this sort is 


Ë Ëp * 
f ar a 3) 
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with g; being the spin degeneracy factor. The Vlasov term 
(the I.h.s. of Eq.(2)) describes the evolution of the distribution 
function in smooth mean field potentials. The collision term 
(the r.h.s. of Eq.(2)) accounts for elastic and inelastic binary 
collisions and resonance decays. The Vlasov term includes 
the effective (Dirac) mass m; = mj +S ;, where Sj = gojo 
is a scalar field; the field tensor F? Y= ôV = avi , where 


Vii = Buju! + Bit p™ + qjA" is a vector field, T° = +1 for p 


and Ā, T° = —1 for p and n; and the kinetic four-momentum 

pe = př- vi satisfying the effective mass shell condition 
*Up* — *2 

pp, =m. 


In the present calculations, the nucleon-meson coupling 
constants g8oN, wN, 8on and the self-interaction parameters 
of the o-field have been adopted from a non-linear Walecka 
model in the NL3 parameterization [23]. The latter gives the 
compressibility coefficient K = 271.76 MeV and the nucleon 
effective mass mą = 0.60my at p = po. The antinucleon- 
meson coupling constants have been determined as 


Sun = —€8un> SpN = E8pn> So = EBON > (4) 


where 0 < é < 1 is a scaling factor. The choice € = 1 corre- 
sponds to the G-parity transformed nuclear potential. In this 
case, however, the Schrödinger equivalent potential becomes 
unphysically deep, Up = —660 MeV. The empirical choice 
of é will be discussed in the following section. 


He = Sea yo a S 
g= Sp + V+ — 
N N Ñ 2my 

The GiBUU collision term? includes the following chan- 
nels: (notations: B — nonstrange baryon, R — nonstrange bary- 
on resonance, Y — hyperon with § = —1, M — nonstrange 
meson): 


e Baryon-baryon collisions: 
elastic (EL) and charge-exchange (CEX) scattering 
BB —> BB; s-wave pion production/absorption* 
NN e NNa; NN e AA; NN © NR; 
N(4, N*)N(4, N*) > N(A)YK; YN > YN; EN > AA; 
EN > AX; EN > EN. 
For invariant energies ys > 2.6 GeV the inelastic pro- 
duction B,B, — B3B4 (+ mesons) is simulated via the 
PYTHIA model. 


Antibaryon-baryon collisions: 

annihilation BB — mesons; EL and CEX scatter- 
ing BB > BB; NN e NA (+ c.c); NN > AA; 
N(A)N(A) > AX (+ c.c.); (DNI) > EE. 

For invariant energies ys > 2.4 GeV (i.e. pip > 1.9 
GeV/c for NN) the inelastic production BB) > B3B, 
(+ mesons) is simulated via the FRITIOF model. 


3 The GiBUU code is constantly developing. Thus the actual version may 
include more channels. This description approximately corresponds to the 
release 1.4.0. 

4 Implemented in a non-RMF mode only. 

5 Described with a help of the statistical annihilation model [24, 25]. 
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e Meson-baryon collisions: 

MN & R (baryon resonance excitations and decays, 
e.g., tN e A and KN © Y*); 2(p)4 e R; AN > aN; 
aN > anN; aN > 74; nN > oN; aN —> ON; tN > 
wraN; nN > rN; n(n,p,w)N —> YK; nN > KKN, 
aN > YKr, nA > YK; KN — KN (EL, CEX); 
KN — KN (EL, CEX); KN © Ya; KN © Y*n; 
KN > EK. 

At Ys > 2.2 GeV the inelastic meson-baryon collisions 
are simulated via PYTHIA. 


Meson-meson collisions: 

M,M, « M3 (meson resonance excitations and 
decays, e.g., nn © p and Kr © K*); MıMı & KK, 
MıMə e KK (+ c.c.). 


HI. ANTIPROTON ABSORPTION AND ANNIHILATION 
ON NUCLEI 


Without a mean field acting on an antiproton, the GiBUU 
model is expected to reproduce a simple Glauber model result 
for the p-absorption cross section on a nucleus (left, Fig. 1): 


> =T tot f dzp(b,z) 
gaber f efi -e > | f (6) 


where Cto is the isospin-averaged total pN cross section. The 


p 


I 
o 
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Fig. 1. Left panel — the straight-line propagation of an antiproton 
in the absence of a mean field. Right panel — an illustration of the 
curved trajectory of an antiproton due to an attractive mean field. 


attractive mean field bends the p trajectory to the nucleus 
(right, Fig. 1). Thus, the absorption cross section should 
increase. 

Figure 2 shows the GiBUU calculations of antiproton ab- 
sorption cross sections on !*C, ? A1 and Cu in comparison 
with experimental data [26-29] and with the Glauber formula 
(6). Indeed, GiBUU calculations without mesonic compo- 
nents of the p mean field, i.e., with scaling factor € = 0, are 
very close to Eq.(6) at Pia» > 0.3 GeV/c. At a lower Plab, 
the Coulomb potential makes the difference between GiBUU 
(£ = 0) and Glauber results. Including the mesonic compo- 
nents of the p mean field, (€ > 0) noticeably increases the 
absorption cross section at Diab < 3 GeV/c. The best fit of 
the KEK data [26] at piap = 470 — 880 MeV/c is reached with 
é = 0.21 +0.03. This produces the real part of the antiproton- 
nucleus optical potential ReV,,, = Us ~ —(150 + 30) MeV at 
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Fig. 2. (Color online) Antiproton absorption cross section on the 
12C, 7 Al, and Cu nuclei vs the beam momentum. The GiBUU 
results are shown by the lines marked with the value of a scaling 
factor £. Thin solid lines represent the Glauber model calculation, 
Eq.(6). For the p+!°C system, a calculation with € = 0 without 
annihilation is shown by the dotted line. 


p = po. The corresponding imaginary part is 


ImVopt = -5 <VpnOtot > P. (7) 

At p = po this gives ImVəp ~ —(100 — 110) MeV indepen- 
dent on the choice of £. It is interesting that the BNL [27] and 
Serpukhov [28] data at pay = 1.6 — 20 GeV/c favor € = 1, i.e. 
ReVop ~ —660 MeV at p = po. This discrepancy needs to be 
clarified, which could be possibly done at FAIR. 

Figure 3 displays the calculated momentum spectra of pos- 
itive pions and protons for antiproton interactions at plab = 
608 MeV/c with the carbon and uranium targets. GiBUU 
reproduces a quite complicated shape of the pion spectra 
which appears due to the underlying 7N < A dynamics. The 
absolute normalization of the spectra is weakly sensitive to 
the p mean field. The best agreement is reached for € = 0.3, 
i.e., for ReVop ~ —(220 + 70) MeV. 


IV. SELFCONSISTENCY EFFECTS 


The strong attraction of an antiproton to the nucleus has to 
influence on the nucleus itself. This back coupling effect can 
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Fig. 3. (Color online) Momentum differential cross sections of 2* 
and p production in p annihilation at 680 MeV/c on °C and 78U. 
The different lines are denoted by the value of a scaling factor £. 
The data points are from [30]. 


be taken into account by including the antinucleon contribu- 
tions in the source terms of the Lagrange equations for w-, p-, 
and o-fields 


O + m0") = Y, TONGO — 8) 


j=N.N 
(O,0" + mp” (x) = >. Eoi TOY PY), (9) 
j=N.N 
dU = 
ð "o (x) + Aa z = 2 PAO: (10) 


with U(c) = imo? + tgo? + 7930", or, in other words, 
by treating the meson fields selfconsistently. As follows from 
Eqs. (4) and (8)-(10), nucleons and antinucleons contribute 
with the opposite sign to the source terms of the vector fields 
w and p, and with the same sign — to the source term of the 
scalar field o. Hence, repulsion is reduced and attraction is 
enhanced in the presence of an antiproton in the nucleus. 
Figure 4 shows the density profiles of nucleons and an 
antiproton at different times when the p implanted at t = 0 in 
the center of the *°Ca nucleus. As a consequence of the pure 
Vlasov dynamics of the coupled antiproton-nucleus system 
(annihilation is turned off), both the nucleon and the antipro- 
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ton densities grow quite fast. Att ~ 10 fm/c the compressed 
state is already formed, and the system starts to oscillate 
around the new equilibrium density pọ = 20. 
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Fig. 4. (Color online) The density of nucleons (thick lines) and 
antiproton (thin lines) as a function of coordinate on z-axis drawn 
through the nuclear center (z = 0). 


Figure 5 displays the time evolution of the central nucleon 
density. The p annihilation is simulated at the time tann. The 
tann = 0 corresponds to the usual annihilation of a stopped p 
in the nuclear center. In this case, the nucleon density remains 
close to the ground state density. However, if the annihilation 
is simulated in a compressed configuration (tanm > 0), then 
the residual nuclear system expands. Eventually the system 
reaches the low-density spinodal region (9 <~ 0.6po), where 
the sound velocity squared c? = 0P/ OP\s-const becomes nega- 
tive. This should result in the breakup of the residual nuclear 
system into fragments. 

A possible observable signal of the p annihilation in a com- 
pressed nuclear configuration is the total invariant mass Miny 
of emitted mesons 


2 
M2, = by pi . (11) 


i 


For the annihilation of a stopped antiproton on a proton 
at rest in a vacuum, Miny = 2my. In a nuclear medium, the 
proton and antiproton vector fields largely cancel each other’. 
Therefore, it is expected that in nuclear medium the peak will 


6 Here, P is the pressure and s is the entropy per nucleon. 
7 The cancellation is exact for the antiproton vector fields obtained by the 
G-parity transformation from the respective proton vector fields, i.e. when 


é=1. 
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Fig. 5. (Color online) The central nucleon density as a function of 
time. The annihilation of p with the closest nucleon into mesons 
is simulated at the time moment tann as indicated. The calculations 
without annihilation and for the ground state nucleus (without p) are 
also shown. 


appear at Min ~ 2m, This simple picture is illustrated by 
GiBUU calculations in Fig. 6. In calculations where tann = 0, 
we clearly see a sharp medium-modified peak shifted down- 
wards by ~ 200 MeV from 2my. The final state interactions 
of mesons create a broad maximum at Miny =~ 1 GeV. For 
annihilation in compressed configurations (fan, = 10 and 60 
fm/c), the total spectrum further shifts by about 100 MeV to 
the smaller Miny. This effect becomes stronger with the de- 
creasing mass of the target nucleus (e.g., for '°O the spectrum 
shift is nearly 500 MeV [14)]). 
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Fig. 6. (Color online) Annihilation event spectrum on the total in- 
variant mass (11) of emitted mesons. Calculations are done for three 
different values of annihilation time tann. 
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V. STRANGENESS PRODUCTION 


Originally, the main motivation of experiments on 
strangeness production in antiproton-nucleus collisions was 
to find signs of unusual phenomena, in-particular, a multinu- 
cleon annihilation and/or a quark-gluon plasma (QGP) forma- 
tion. In Ref. [31], the cold QGP formation has been suggest- 
ed to explain the unusually large ratio A/K} ~ 2.4 measured 
in the reaction p!*!Ta at 4 GeV/c [32]. On the other hand, 
in Refs. [8, 11, 16-18, 33-35] most features of strangeness 
production in pA reactions have been explained by hadronic 
mechanisms. 

Figure 7 presents the rapidity spectrum of (A + 5°) 
hyperons, K? mesons and (A + È?) antihyperons for the col- 
lisions p(4 GeV/c)!*!Ta in comparison with the data [36] and 
the intranuclear cascade (INC) calculations [11]. The GiBUU 
model underpredicts hyperon yields at small forward rapidi- 
ties y ~ 0.5 and overpredicts K? yields. In the GiBUU calcu- 
lation without hyperon-nucleon scattering, the (A + °) spec- 
trum is shifted to forward rapidities. However, the problem 
of underpredicted total (A + X?) yield remains. A more de- 
tailed analysis [16] shows that 72% of Y and Y* production 
rates in GiBUU are due to the antikaon absorption processes 
KB — YX, KB — Y*, and KB — Y*n. The second largest 
contribution, 23% of the rate, is caused by the nonstrange me- 
son - baryon collisions. The antibaryon-baryon (including the 
direct pN channel) and baryon-baryon collisions contribute 
only 3% and 2%, respectively, to the same rate. The under- 
prediction of the hyperon yield in GiBUU could be due to the 
used partial KN cross sections, in-particular, due to the prob- 
lematic K~n channel®. The possible in-medium enhancement 
of the hyperon production in antikaon-baryon collisions is al- 
so not excluded. 

As shown in Fig. 8, at higher beam momenta the agree- 
ment between the calculations and the data on neutral strange 
particle production becomes visibly better. The exception is 
again the region of small forward rapidities y ~ 0.5 where 
both GiBUU and INC calculations underpredict the (A + % 0) 
yield. 

Finally, let us discuss the £ (S = —2) hyperon production. 
The direct production of Æ in the collision of nonstrange par- 
ticles would require to produce two s5 pairs simultaneous- 
ly. Thus, £ production could be even stronger enhanced in a 
QGP as compared to the enhancement for the § = —1 hyper- 
ons. Fig. 9 shows the rapidity spectra of the different strange 
particles in p!” Au collisions at 15 GeV/c. Even at such a 
high beam momentum, the S$ = —1 hyperon spectra still have 
a flat maximum at y ~ O due to exothermic strangeness ex- 
change reactions KN — Yr with slow K. In contrast, the 
second largest, ~ 18%, contribution to the production is 
given by endothermic double strangeness exchange reactions 
KN — EK °. Since the threshold beam momentum of K for 


8 The K~n channel has been improved in recent GiBUU releases, however, 
after the present calculations were already done. 
9 The main, ~ 24%, contribution to the total yield of 5’s at 15 GeV/c is 
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Fig. 7. (Color online) Rapidity spectra of (A + X°), K®, and (A +5”) 
from p'*'Ta collisions at 4 GeV/c. See text for details. 


the process KN —> &K is 1.05 GeV/c, which corresponds to 
the KN c.m. rapidity of 0.55, the rapidity spectra of &’s are 
shifted forward with respect to the A rapidity spectra. How- 
ever, in the QGP fireball scenario [31], the rapidity spectra of 
all strange particles would be peaked at the same rapidity. 


VI. SUMMARY 


This work was focused on the dynamics of a coupled 
antiproton-nucleus system and the strangeness production in 


given by Æ* — Ær decays. The direct channel NN — 5£ contributes 
~ 10% only. 
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Fig. 8. (Color online) Rapidity spectra of A, K? and A from p*Cu 
collisions at 8.8 GeV/c. The data and INC calculations are from 
Ref. [33]. 


pA interactions. The calculations were based on the GiBUU 
transport model. The main results can be summarized as: 


e The reproduction of experimental data on pA absorp- 
tion cross sections at Plab < 1 GeV/c and on x* and p 
production at piap = 608 MeV/c requires a strongly at- 
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tractive pA optical potential, ReVəp =~ —(150 — 200) 
MeV at p = po. 


e As a response of a nucleus to the presence of an an- 
tiproton, the nucleon density can be increased up to 
p ~ (2 — 3)po locally near p. Annihilation of the p in 
such a compressed configuration can manifest itself in 
the multifragment breakup of the residual nuclear sys- 
tem and in the substantial (~ 300—500 MeV) shift of an 
annihilation event spectrum on the total invariant mass 
of produced mesons Miny toward low Miny. 


e GiBUU describes the data on inclusive pion and proton 
production fairly well. Still, the strangeness produc- 
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Fig. 9. (Color online) The rapidity spectra of (A + 2°), K?, (A + °), 
Æ, and 2° from p!®’Au collisions at 15 GeV/c. 


tion remains to be better understood (overestimated K? 
- and underestimated (A + 2°) - production). 


e = hyperon forward rapidity shift with respect to A is 
suggested as a test of hadronic and QGP mechanisms 
of strangeness production in pA reactions. 
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